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Abstract 

The formalism based on the equal-time Wigner function of the 
two-point correlation function for a quantized Klein-Gordon field is 
presented. The notion of the gauge-invariant Wigner transform is 
introduced and equations for the corresponding phase-space calculus 
are formulated. The equations of motion governing the Wigner func- 
tion of the Klein-Gordon field are derived. It is shown that they lead 
to a relativistic transport equation with electric and magnetic forces 
and quantum corrections. The governing equations are much simpler 
than in the fermionic case which has been treated earlier. In addition 
the newly developed formalism is applied towards the description of 
spontaneous symmetry breakdown. 
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1 Introduction 



We consider the time evolution of quantum fields interacting with a classical 
(external) field. A useful example of such a system is the particle creation 
by strong but slowly varying electromagnetic fields. In this case the particle 
fields require full quantum treatment while the electromagnetic field may be 
considered classically. It is sometimes necessary to treat these classical fields 
as fully dynamical objects rather than as a fixed background. The study of 
the back reaction problem is such Quite recently this type of dynam- 

ical description has been successfully applied towards the description of 
the electron-positron vacuum in QED. Eventually some of these approachs 
may serve as crude models for the formation stage of the quark-gluon plasma. 

To achieve their goals the authors of had to simplify the notion of the 
vacuum by truncating the infinite hierarchy of correlation functions. While 
the full description of a vacuum state at a given instant of time is pro- 
vided by the set of all possible n-point equal-time correlation functions of 
the quantum field, two-point functions are sufficient to construct the electric 
current and therefore to describe the influence of the electron-positron vac- 
uum on a classical electromagnetic field. The authors of [0J make use of the 
Fourier-transformed two-point equal-time correlation function (an analog of 
the Wigner function) for fermionic fields. The resulting equations describe 
the evolution of the vacuum from a phase-space point of view. Our goal 
here is to extend the basic results presented in the cited work to the case of 
scalar electrodynamics. The relative simplicity of the scalar result helps in 
the interpretation of the physical features of the interacting system. 

The use of the Wigner functions to the description of the QFT vacuum is 
not new P, [|, H. Neither is the use of the two-point correlation function to 
describe the interaction with a classical but dynamical electromagnetic field 
[@- 0' HI • The main difference is that we use equal-time correlation functions 
and do not perform the Fourier transform with respect to the time variable. 
(The equal-time approach is also used in Ref. 0.) The equation of motion 
is therefore a true evolution equation that can be numerically solved as an 
initial-value problem. 

In section |2] we introduce the gauge-invariant Wigner transform and show 
how the noncommutative operator algebra of quantum mechanics is imple- 
mented in the space of phase-space transforms. The Wigner transform of 
the equal-time two-point correlation function of the Klein-Gordon field in 
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Feshbach-Villars representation is defined and discussed in section [3], which 
will, incidentally, illustrate the power of the methods developed in section 
|^. Finally, in section [| we apply this formalism to spontaneous symmetry 
breakdown. 

2 Wigner function calculus 

The Wigner transform associates a quasiprobability distribution in classical 
phase space to every quantum-mechanical density matrix in such a way that 
quantum-mechanical averages can be computed like classical statistical av- 
erages. The equation of motion of the quasiprobability distribution derives 
from the quantum-mechanical Hamiltonian and corresponds to the classical 
Liouville equation with quantum corrections which can be given as an ex- 
pansion in h. Since the Liouville equations is easier to solve than the full 
quantum treatment, the Wigner function provides a tool to simplify calcula- 
tions when semiclassical approximations are valid. 

Similarly, the two-point function of a quantum field theory can be inter- 
preted as the matrix element of a density matrix. The resulting phase-space 
quasiprobability distributions reflect the many-particle content of the quan- 
tum field theory. Their equation of motion follows analogously from the 
operator equation of motion, i.e., the wave equation, and the corresponding 
Hamiltonian. 

Mathematically, the Wigner transform provides a realization of the ab- 
stract Heisenberg operator algebra of quantum mechanics on the space of 
functions over phase space, thus translating every operator equation — like 
the Heisenberg equation of motion — into a differential-equation on phase 
space — like the Liouville equation. Especially, it induces a noncommutative 
multiplication on this function space that is characteristic of the quantum- 
mechanical content of the theory. Once the representation of the basic op- 
erators and this multiplication is known, every operator equation can easily 
and without explicit calculation be translated into a phase-space equation. 

In the following, we discuss this calculus for a Wigner function exhibiting 
local U(l)-gauge invariance. The notational conventions are from the review 
article fllQ| , where the gauge issue was not considered. We set c = 1 but keep 
h in this section since expansions around the classical limit are frequent. 
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2.1 Definition and elementary properties 

Let us consider a D- dimensional configuration space of vectors q G R/° and 
the corresponding Hilbert space Tt, spanned by basis vectors (q\. Suppose 
that the operator A transforms covariantly under a local U(l)-gauge trans- 
formation A (9) 

(q\A\q') -> e ie{A{q) - A{q ' )} (q\A\q'). (1) 
The corresponding transformation of the vector potential A(q) G is 

A(q)^A(q) + ^(q). (2) 

Every operator A may be associated with a function in phase-space (q, p) G 
R°xR° 



/(% %€. ri^vl^ \ 

d °y (q - \y\ A \q + \y) exp f -p • y + — J _ ^ A(x) ■ dx 1 . 

" (3) 

This is an obvious generalization of the well-known Wigner transformation, 
and we will use this name here. Indeed, Pp(q,p) (p being a density matrix) 
is the usual Wigner phase-space distribution. The line integral in Eq. is 
taken over a straight line, i.e., 

pq+y/2 p+1/2 

/ A(x) ■ dx = / dXy ■ A(q + Xy) (4) 

Ja-v/2 J -1/2 



where the dot indicates the scalar product in D dimensions. This factor 
makes the resulting function P^(q,p) a gauge-invariant realization of the 
abstract operator A in phase space. Every operator equation can thus be 
translated into a phase-space equation. In the following we discuss the re- 
quired formalism and apply it to the Heisenberg equation of motion for a 
density matrix. 

The inverse of the transformation (^) is 

f d^ k f % zc fi+y 

(q-y\A\q + y) = J J^yjPA^ k) exp (-2-* • y - - A(x) • dx 

(5) 

The matrix representation and the phase-space representation of the operator 
A are therefore equivalent. 
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The Wigner transform ([3]) conserves the linear structure of the underlying 
operator algebra 



Pa+b(q,p) = Pa(q>p) + p&(q,p), p*a = ®Pa- (6) 

If the operators carry internal degrees of freedom, the Wigner function will 
acquire a matrix structure, and the second equation in (|j) holds for con- 
stant matrices a. Because of noncommutativity, multiplication of operators 
translates into a more complex phase-space operation considered below. 

The Wigner transforms of the position operator q and the momentum 
operator ft can be determined directly from their matrix elements, 

d 

(q\q\q') = 5{q - q), (q\p\q'} = -ih8(q - q) — - - . (7) 

For the position operator, insertion into (||) yields 

P<i(q,p) = q- (8) 

This relation may be generalized to analytic functions W(q): Pw(q)(q,p) = 
W(q) by a Taylor expansion. In the case of the momentum operator, the 
derivative in (|7]) gives an additional term from its action on the line integral, 
which yields 

Pp(q,p) =P + eA(q). (9) 
Therefore, the kinetic momentum 7r = p — eA{q) has the Wigner transform 

P*(q,p)=p. (10) 

This allows us to interpret the coordinate p as a true (i.e., kinetic) momen- 
tum. 



2.2 Noncommutative multiplication in phase space 

The noncommutative multiplication of the operator algebra translates into a 
'noncommutative multiplication' of their Wigner transforms. Consider 



Pab^p) = J d D y(q~ \y\AB\q + \y) 



x exp \ —p ■ y + — / A(x) ■ ax . (11) 

\h k Jq-y/2 J 



Inserting a complete set of intermediate states and using (|5|) gives us — after 
some shifts of integration variables — the multiplication formula in its integral 
form: 

Pab&p) = W 2 °J d D yd D y'd D kd D k' 

x Pa(q + y> p + k ') p B (<? + y',p + k) 
x exp {-jr{yk- iA')} 

{ie ( ri-y+v 1 rq+y-y' rq+y+y' \ 1 
T / + / + / A(x) -dx\. (12) 

fi \ Jq+y-y' Jq+y+y' Jq-y+y' ) J 

With A(x) put to zero one gets the multiplication formula given by JTO ] . 

Let us now consider the usual case D = 3 with A(x) being a three- 
dimensional vector potential. The closed-line integral in (|I^ )is related to 
the surface integral over the magnetic field B(q): 

( A(x)-dx= f d 2 xn-8(x), B(q) = curl A(q), (13) 

JdA J A 

where A represents the integration region (here a triangle) and dA its border, 
n is the normal to the triangle plane, dx the line element of the border, and 
d 2 x the area element of the triangle. We parametrize the surface of the 
triangle according to 

x = q + Aiy + X 2 y', -1<Ai,A 2 <1, Ai + A 2 > 0. (14) 
The surface integral in question may be rewritten: 

/ d 2 xn -cm\A{x) = C dAi C d\ 2 (y x y') ■ B{q + X x y + \ 2 y'). (15) 
Ja J-i J-\i 

Inserting the above formulae into ([12]) we get 

P AM (q,p) = (^)- 6 | d 3 yd 3 y'd 3 kd 3 k'P A (q + y,p + k') 

( 2i ic f 

x exp |-(j/fc - y'k') + jr{yxy')- J dA x dA 2 i3(g + \ x y + \ 2 y') 



( , d ,0 s 
x exp ?/ • — + k ■ — 
\ oq op 



(16) 

q=q,P=p 
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where we made use of the 'translation' formula 

P B (q + y',P + k') = exp (y' ~ + k ■ P 6 (q, p) 



(17) 



q=q,p=p 



Now the integration over k can be performed, yielding a (formal) delta 
function 

Subsequent integration over dy results in the expression: 

exp f ee -4|^ L dXidx * B ( q - A 4| + xw ) } 



(19) 



q=q,p=p 

We can now get rid of the y' in the first exponential by substituting it by 

% d 
Yidk r 



y 



(20) 



where the derivative acts on the last exponential. This leads to the expression 
4, ^d Pl dk> L dXldX2Bk \ q - Xl 2ldp- X Ydk')) 



(eh d d 
x exp < ~—e ijk - 



(21) 



q=q,p=p 

After partial integration, this operator acts on P^ and can be replaced by 

—d /dp, where the arrow indicates that the momentum derivative acts on 
the function to its left yielding the final expression 

J., / s t> ( h 9 h d \ 



eft d d I h d h d 

x exp { -—tijk — — J i dAi J d\ 2 B k I q - Ai — — + A 2 



4z (9pj j-i 7-Ai I 2i dp 2i dp 

In this way the left multiplication with an operator A was translated into a 
phase-space operation. For the right multiplication we get 

{eh d d r 1 /— A i „ „ / . ft 9 , h d 
Tre^-fc— — / dAi / d\ 2 B k q - Ai — — + A 2 — ^~ 
4« dpi opj J -i J -i \ 2i op 2i op 

x^fo#U,p= p - (23) 
2.3 Kinetic energy and the flow term 

In the following we will derive equations of motion for the Wigner transform 
of the density matrix p. The time evolution of the latter is governed by the 
Heisenberg equation of motion 

i = 4[*.fl- (*) 

Therefore we have to find the phase-space analog of [H, p] . We do not as- 
sume any special properties of p, and the following would apply to a general 
Heisenberg operator that has no explicit time dependence. (Allowing for an 
explict time dependence does not pose any problem.) 

We begin with the kinetic momentum ft = p—eA(q). Its Wigner transform 
is Pjt{q,p) = P- In this case only first-order derivatives d /dp in (^) and ( p3|) 
lead to a nonzero result. Therefore only the first two terms in the Taylor 
expansion of the exponential are to be left 

!eh d d /"i / Ti d HO 
— e Hk-^ — t; — / dAi / dA 2 £>fc [q + Ai — — A 2 — 77- 
Ai dpi dpj J-i j-Ai \ 2% dp 2i dp 

- 1+ f £ '4£/- 1 l dAl /l dAA ( 9+Al s|) < 25 > 
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are to be left. Employing this we find 

h d ieh r 1 

2idq 4 J-i J-x 1 " " z ~ \^ 2i dp J dp 

P P (q,p), (26) 

(27) 

From this one easily obtains the following formulae for the commutator and 
anticommutator of n and p: 

\ d ieh f 1 , , M / h\ d \ d ) „ . w . 

^fori = {-*aj-^-y. I d «(»-ar*j x *} p ' ( '' 1 '» 28) 
WwO = {*-^/ 1 <^(«-£|)*|}'V«.i>>- < 29 » 

It is convenient to define two phase-space operators /) and P : 

P> = — + e / dAB U + ^A— x — 
[ Oq J -i/2 \ op J op) 

" T q +eB(q) *h (30) 

dAAS g + i/iA— x — 
-i/2 \ apy 

where the expansion to lowest order in H is given. These operators act on 
the Wigner transforms in a nonlocal manner. They are generalizations of the 
g-derivative and the p-multiplication in phase-space and their use simplifies 
our notation: 

P[*A<1,P) = ~ihDPp(q,p), P{*,p}(q,p) = 2PP p (q,p). (32) 

We illustrate his technique using a simple example. Let us consider the 
Hamiltonian of a nonrelativistic particle coupled to an external electromag- 
netic field 

H = H kin + eA°(q), (33) 

^ _ tf^.^^.i, (34) 
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|p - ieh J 1 ' dXXB [ q + ihXj^ ] x £- \ w p, (31) 



where A°(q) is the time-like part of the electromagnetic four-potential. 

To study the time evolution of the Wigner function P p (q,p), we must 
evaluate the commutator and, in the case of a more complex Hamiltonian 
with matrix structure, the anticommutator. Using the relations found above 
this is straightforward. We rewrite both as 

[n 2 ,p] = {*,[#, p)}, (35) 
{vr 2 ,/)} = fop}} + (36) 

The double (anti) commutators correspond to a repeated application of the 
operators (0) and (|3l| ): 



P ■ D 

P\^/2m,p]{<i,v) = -ih Pp(q,p), (37) 

P{nV2m,p}(^P) = ^[p 2 -\D 2 ^P p (q,p). (38) 

In the expression for the commutator we note the classical Liouville flow term 
with the minimal coupling to a magnetic field, 

P-D P f 9 „, , d\ 
m m \oq op J 

where only the lowest order term in fi was left. 

The commutator involving a scalar potential U(q) can also be rewritten 
clS cL differential operator in phase space 

rl/2 f)TT ( f) \ f)TT 

r [ vd«>P)= i J_ 1/2 ^{* + ih %) p ^ ( 4 °) 

The correspondig expression for the anticommutator is 

« C/(?)Pp(?,p)+0(S 2 ). (41) 

The time evolution of a Wigner function is not only determined by the 
change of the underlying operator, which is usually given by the Heisenberg 
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equation of motion, but also by the change in the line integral due to a 
time-dependent vector potential. The additional term can be rewritten in a 
manner similar to that used above: 

!%p) = JV, + «££m(* + «a|;) -|p,. (42) 

Inserting the Heisenberg equation of motion (|24|) with the Hamilton operator 
of Eq. (^), we obtain 

D t P ii (q,p) = -!—^Pe(q,p) (43) 
m 

with the generalized time derivative 

d r 1 / 2 ( d \ d d d 

D '-m + U eS {" + m T P ) ■ ai *> at + e£(9) ' 55 + ' ' ' (44) 

The derivatives of A resulting from ( |i0|) and (^) were combined to yield an 
electric field 

£ ^ = - a 4 - f • < 45 » 

To lowest order in this is the nonrelativistic Liouville equation of a charged 
particle in an external electromagnetic field. 



2.4 Calculation of observables 

The Wigner function is normalized to give 

J dD Vj^ P A^P) = Jd D p(q\A\q)=TrA. (46) 

In the case of a density matrix it is just one. Observables are computed by 
taking the trace over the product of density matrix and observable operator: 
Using the multiplication formula ([12]) one can prove that 

(O) =Tr (Op) = J d D q^ B TT(P 6 (q,p)P i) (q,p)), (47) 

where Tr is the trace over the matrix structure of P(q,p) (if necessary). This 
formula allows us to compute any observable from the phase-space repre- 
sentation of the density matrix; it further suggests the identification of the 
phase-space function of an observable operator and the classical observable 
function. 
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3 The Wigner function of the Klein— Gordon 
field 



We will now treat the equal-time Wigner function of a Klein-Gordon field 
in an external electromagnetic field. Instead of the matrix elements of a 
density matrix, the object to consider will be the symmetrized two-point 
correlation function. Its properties are very similar to a quantum-statistical 
density matrix in that it describes a many-particle system. All single-particle 
observables can be computed from the two-point function, and thus from its 
Wigner transform. 

The time evolution of the two-point function is derived from the Klein- 
Gordon equation for the Heisenberg field operators. Similarly to the nonrel- 
ativistic problem treated above the Klein-Gordon equation can be translated 
into a differential equation in phase-space with its leading contribution being 
just the classical relativistic flow equation. Two complications arise: The 
Klein-Gordon equation contains a second-order derivative in time. Its trans- 
lation into phase-space would also contain a second-order derivative which 
does not lead to a sensible evolution equation. The reason for this is that 
the Klein-Gordon field has an internal charge degree of freedom which can 
be made explicit in the two-component formalism of Feshbach and Villars 
[[LlH . In this representation, the time evolution is first-order in time and gov- 
erned by a nonhermitian Hamilton operator with a 2 x 2-matrix structure. 
As a second complication in a fully interacting case (scalar electrodynamics) 
the equation of motion of the two-point function involves higher correlation 
functions. Here we consider the simpler problem of the interaction with a 
classical albeit dynamical electromagnetic field in which case the equation of 
motion for the correlation function is closed and the electromagnetic field is 
governed by the Maxwell equations. 

3.1 Definition 

The Klein-Gordon field which obeys the wave equation (h = c = 1) 

((<9 M - ieA^) {&* - ieA^) + m 2 ) <j>{x) = 0, (48) 
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can be expressed by a two- component wave function 




) 



(49) 



transforming regularly under U(l)-gauge transformations where 




m dt m 



i dip eA l 



) 




if, % d<p eA { 



(50) 



The equation of motion for the two-component field operator is: 




in the first-order giving rise to a Schrodinger type evolution equation. The 
right-hand side can then be interpreted as a Hamiltonian operator acting 
on <3> and the formalism of the previous section can be utilized to derive a 
phase-space equation. The fact that the resulting Hamiltonian operator is 
not hermitian does not cause any severe problems, since the nonhermiticity 
is confined to the matrix structure. 

In the following, we consider the time evolution of the symmetrized cor- 
relation function 



which is a 2 x 2 matrix. The matrix element is taken in the Hilbert space of 
quantum field theory with respect to some state Q. Later on this state will 
be identified with the vacuum. We may consider C + (q,q') to be the matrix 
element of some density matrix p in the abstract Hilbert space spanned by 
position eigenvectors (q\, 



CUq,q';t) = (n\{$ a (q,t),$+(q',t)}\n) 



(52) 




(53) 



By Eq. (|5"T|), the density matrix p evolves under the Heisenberg equation of 
motion 



i^ = Hf>- pH + (54) 
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with the two- component nonhermitian Hamiltonian 



H 



(p - eAf 
2m 



1 1 
-1 -1 



+ m 



1 
-1 



+ eA°l 



(55) 



where p denotes momentum operator —id/dq. 

The energy, momentum, charge, and current of the Klein-Gordon field 
are bilinear forms in the field operators that can be computed from the two- 
point correlation function or from the Wigner function. For example, the 
energy and momentum of free particles are given by the elements of the 
energy-momentum tensor, 
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oo 



2m \ dt 



dt 



e 



Oi 



1 1 f 

1 1 ) 2m 



(V0*)(V0) +m 2 <j)*<f) 



1 
1 



m 



-m 



1 1 
-1 -1 



V 



(56) 
(57) 
(58) 
(59) 



The sum of these quantities and the corresponding quantities of the electro- 
magnetic field are conserved. Charge and current density operators read 



J 



le i 
2m \ 

e 
2mi 



dt T dt J 



m 



e 2 A 



m 



e A 

— P 

m 



1 
-1 

f 1 1 

1 1 



$, (60) 
$. (61) 



The expressions in brackets constitute the observable operators in the sense 
of Eq. (M. 



3.2 The Wigner function and its time evolution 

The Wigner transform associates a matrix-valued phase-space function 
P(q, p) with the density matrix p. Its evolution is governed by Eq. ( [54] ) where 
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the nonhermiticity of the Hamiltonian (|55| ) has been taken into account. This 
may be translated into a phase-space evolution equation. Defining 

a=(a 3 + ia 2 ) =(li _\ J ' b = a 3 = ( I 1 (62) 
and using ( |4"2| ) one obtains: 

P Hp-pH+ (?, P) = ( a ' P [J5 2 -P] (?> P) + P [p 2 ,p] (?> P) ' G+ ) 

+ i ( G ' Apv^p) ~ p {P 2 ,p}(?'P) • a+ ) 

+ m(b- P{q, p) - P(q, p) ■ b + ) . (63) 



Following the general procedure outlined in the previous section we get the 
phase-space equation of motion: 

iD t P(q,p) = -i—-b(a-P(q,p) + P(q,p)-a + ) 

+ 1 (f* - Id 2 ) (a ■ P{q,p) - P(q,p) ■ a 

+ m(b- P(q, p) - P(q, p) ■ b) . (64) 

It is convenient to expand P{q,p) (which is a 2 x 2 matrix) over the 
Pauli matrices (7j and the unit matrix, assembling the coefficients into a 
four-component object / , 

P(q,p) = fs(q,p) I + £/ 3 -i(<Z,P,tK (65) 

i=i 

The equations of motion take the form 

P ■ D 

A/o = / 2 + / 3 , 

m 



A/i = -(-^D 2 --P' 2 )(f 2 + f 3 ) + 2mf 2 , 
\4m m J 

A/a = — — W 1 + / -2m/ 1; 



1 1 



4m m / m 



A/a = -(^-&-—)fi-—f . (66) 

V 4m ml m 
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The components of / may be interpreted as the phase-space densities of 
simple observables: 

charge: ef , (67) 
p 2 

energy: — (f 2 + f 3 ) + mf 3 , (68) 
2m 

P 

current: — e(/ 2 + / 3 ), (69) 
m 

momentum: p(f — fx). (70) 

It may be interesting to note the conservation law that follows from the 
preservation of the norm 

9 ' dD?d >(/„ 2 -/?-/l + /3 2 )=0. (71) 



dt J (2irh) 3 

The Feshbach-Villars Hamiltonian (EST) does not give rise to unitary evolution 



but to r-unitary evolution |L2] so that two fs come in with the opposite sign. 
The conservation laws for energy-momentum and currents may be easily 
derived and we do not state them here. 

The free vacuum solution can be derived from the explicit expression for 
the correlation function C + : 

fo = 0, h = 0, 

/ 2 + / 3 = ^, f 3 ~f2 = ^. (72) 
hp m 

It is the stationary solution of (|66|) . The vacuum has a constant energy den- 
sity (which must be subtracted when computing observables) but no charge. 



3.3 Diagonalization of local oscillations 

For a further discussion of these equations we consider the case of a electric 
field slowly varying in space. This will allow us to drop the higher-order 
corrections in P and D. Nevertheless, the equation of motion ( |S6D contains 
self-couplings of the components that make its solutions nonstationary even 
in the absence of a gradient. The resulting oscillations are related to the 
internal charge degree of freedom. The eigenfrequencies of these oscillations 
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are (in the absence of an electromagnetic field) just rk2E p , showing clearly 
that they are related to relativistic effects and therefore to the existence of 
antiparticles. Let us expand the Wigner function P(q,p) into eigenmodes Pj 
of these oscillations: 

P(q,p) = f2fi p i, (73) 

i=i 



with 



Pi = 

Pi = 



m/Ep + E p /m m/ 'E p — E p /m 

m/Ep — Ep/m m/E p + E p /m 

1 0\ 
1 ) ' 



~ if m/Ep - Ep/m m/E p + E p /m =f 1/2 \ ( . 

3/4 " 4 V m / E P - E P / m ± 1/2 - ^ P / m / ' 

The first two matrices are connected with the eigenfrequency zero while the 
last two correspond to ±2E P . Each of these matrices can be constructed 
from plane-wave solutions of the Klein-Gordon equation. Note that the 
basis matrices Pi depend on the energy E p and are therefore different in 
different parts of the phase space. A particle whose inner degree of freedom 
is oscillating in an eigenmode of these equations will be out of phase when it 
gains or loses momentum. 
With the approximation 

1 d 2 

p 2 + m 2 + -— ^p 2 + m 2 , (75) 
4 o 2 q 

which is justified in the semiclassical case, the equations of motion take on 
the form 

*<«>•£) - -i- & ™ 



_/ 2 + ££(9) ._j = _2L._ ( / 1 + / , + /l ), (77) 

d ~ d \ 1 p d ~ 

dt h + e£{q) -dp-) = -2Yp-d- q h 
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E 2 ~2~ p ^ 3 ' ^ ^ 



|/ 4 + e^)-|) = -^-| /2 



+ e ^t±-2iE p U (79) 



Here fi, f 2 and the linear combination / 3 + / 4 are real, while / 3 — / 4 is purely 
imaginary. 

These variables have the meaning of the following phase-space densities: 

charge: ef 2 , (80) 
energy: E p f x (81) 
current: ^(f 1 + f 3 + f i ), (82) 

momentum: — pf 2 + ip{fz — fi)- (83) 

The Wigner transform of a positive- or negative-energy plane-wave solution 
with the momentum p of the Klein-Gordon equation is 

fi = lS(p-po), f 2 = ±15(p-p ), f 3 = h = 0. (84) 

Wigner transforms with nonzero / 3 and / 4 involve interferences between 
positive- and negative-energy modes ("Zitterbewegung"). The correlation 
function of a free vacuum may be expressed as the superposition of all 
positive- and negative-energy contributions (with the same sign, as opposed 
to the Feynman propagator). The free vacuum solution is therefore 

A = l, h = h = h = 0. (85) 

In Eq. (ff6|)— (79|), the classical limit is exhibited in a more elegant way. 
The terms on the right-hand sides of Eqs. ( [T8D and ([79]) are comparable only 
if the variation of f 2 is on a scale of E p , i.e. the Compton wavelength. In 
a classical limit, where this is not the case, the last two equations decouple 
since / 3 and / 4 will vary so fast that their contribution to f 2 is averaged out. 
By introducing two new quantities 

/+ = §(/!+ A), /! = /+ + /-, 

/- = |(/i-/ 2 ), h = /+-/-, (se) 



18 



which obey 



d p d 

di f+ = ~E- p 'dq~ U ' (8?) 

y- " < 88 » 

one obtains two decoupled Vlasov equations of motion. As can be seen from 
(|80"D, /+ and /_ carry different charges and thus describe the phase-space 
densities of positive and negative particles. In this case we have a collisionless 
relativistic gas consisting of positively and negatively charged particles. 

Going back to the equations (1761)— ([79) we can give some meaning to the 
individual terms. As we have seen fi and f% are the particle and charge 
number density and therefore represent the 'substantial' part of the Wigner 
function. Quantities f% and arise from interferences between these two 
parts. On the right-hand side, the local derivatives involving f% and cor- 
respond to flow while those involving and f± implement local interference 
(together with the neglected second derivative). Most important, the terms 
proportional to 

(89) 

stem solely from the momentum dependence of the basis ([71|). They give rise 



to pair creation fl3| . This will be discussed in detail in a future publication 



4 Symmetry breaking 

As mentioned in section [I], the phase-space description may be useful to study 
the back reaction problems. We would like to ilustrate this point with a rel- 
atively simple example. Namely, we now show that the formalism developed 
above provides a useful tool for the dynamical description of spontaneous 
symmetry breaking. The Klein-Gordon equation coupled to a scalar poten- 
tial reads 

{id t ) 2 <j)= j <f) + m 2 <f) + U{q,t)(f). (90) 
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This leads to the interaction Hamiltonian (in Feshbach-Villars form) 

U{q,t) 



H, 



int 



2m 



1 1 
1 1 



Using 



and PT| ) and expanding as in 




we get 



dt 



fo 



-~{f2 + fs) + ^-U.{q) (h + fs), 
m Oq 2m 



9 f 



1 d 2 



I) 



4m dq 2 m ) 



(f 2 + f 3 ) + 2mf 2 



d_ 

dt 



dt 



+ ^U + (q)(f 2 + f3 ), 

4m oq z m J m oq 

---L{U^q)U + U + {q)h), 
2m 



4m dq 2 
1 



p 2 \ p d 

— Ji ' tt/o 

ml m oq 



f 2^( U -(l)fo + U+(q)fi), 



with the potentials U + and LL given by 



2dp' 
2U(q,t)+0(h 2 ), 



i d 



U-(q,t) 



-rr ( i d \ rr ( id 



(91) 



(92) 



(93) 



(94) 



(95) 



(96) 

+ O(h). (97) 



We shall consider the case when the potential U is generated by the 
scalar field itself. To illustrate this point let us consider the charged field (f> 
described by the Klein-Gordon equation (|90|) with 



U 



-M 2 + \\d)\ 2 



(98) 



where the positive constant A defines the strength of |0| 4 self-coupling. 
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It is quite easy to find all time (and coordinate) independent c-number 
solutions to this equation. For M 2 < m 2 the only possibility is 



= 



(99) 



while for M 2 



> m 



2 there is an additional set of solutions 



|0| = \j(M 2 - m?)/\. 



(100) 



The solutions (|100[ ) are known as the symmetry-breaking solutions. 

It is interesting to calculate the Wigner function corresponding to these 
classical solutions. The case = leads to the null Wigner function and is 
of no special interest. On the other hand the symmetry breaking solution 
gives rise to 



while all other components of the Wigner function are equal to zero. It is 
easy to show that the above Wigner function satisfies our basic equations 
fl92"D-( |95[) with the scalar potential given by 



which follows from Eq. fl98D and the definitions of f2 and f^. Note that this 
form of potential introduces self-coupling of a mean-field type. 

The appearance of the Dirac delta function reflects the infinite-range cor- 
relation of the classical solution ( P-OOQ . The infinite-range correlation (of 
quantum fields) indicates symmetry breaking. 

The classical considerations presented above do not allow for the dynam- 
ical description of a symmetry-breaking phase transition. Let us consider a 
simple — albeit rather abstract — situation in which one is allowed to switch 
on the — M 2 in fl9~8| ) at a given time t = t - (We will assume that it was 
switched on instantly but other possibilities are not excluded.) Initially the 
quantum field is in its symmetric vacuum state with mass equal to m. 
After the external potential — M 2 is switched on this will no longer be a 
stable state (for suffiuciently large M. Classically one has to provide a small 
and quite arbitrary deviation from = to initialize the fall towards the 
"true" vacuum. There is also no way to damp the resulting oscillations. 




(101) 




(102) 
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Therefore it is not possible to predict the average time necessary to complete 
the phase transition from purely classical considerations. On the other hand 
quantum theory provides the fluctuations necessary to drive the transition. 
Indeed, equations d92|)- (|95[) contain the necessary quantum fluctuations and 
it is possible to study the dynamics of symmetry breaking using them. 

For simplicity we concentrate on the spatially uniform and isotropic solu- 
tions. In that case functions / depends only on \p\; they do depend neither 
on q nor on the direction of the momentum p. The potential U is also uniform 
and our evolution equations (|9^)-(|95|) reduce to: 



§- f fo = 0, (103) 

= -(/ 2 + /s) + 2m/ a + ^(/ 2 + /s) ) (104) 

at m m 

777 /2 = fx - 2m/i /1, 105 

at m m 

777/3 = — h H h- (106) 

at mm 

Note that f decouples completely from the rest of the system and does not 
undergo any evolution. 

The time dependent scalar potential is chosen to be: 

U(t) = -M 2 9(t -t o )+\J-0-(f 2 + f 3 ) + C. (107) 

The first term is some external potential that is switched on (rapidly) at 
t = to- If M 2 > m 2 it will change the sign of a "mass squared term" inducing 
the transition towards new equilibrium. The second term describes a mean 
field self-coupling of |</> | type. It is well known that for A > this type of 
interaction will tend to stabilise the system. In addition we have a constant 
C; its value will be specified in a moment. 

As an initial condition we assume the free vacuum state with distribution 
functions given by Eq. (f72|) i.e., 



/o = -°° = 0, ft' 00 = 0, 

ft= — 00 1 ft= — 00 ^ ft= — oo rt= — oo V /iao\ 

h + h = TT> h ~ h = — » l luy J 
hip m 
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were superscript t = — oo tells us that we assumed this form of Wigner 
function in a distant past. Note that due to the self-coupling term this state 
is no longer stable (in general) and will undergo some evolution even for 
t < to i.e., in absence of the external potential —M 2 Q(t — to). This means 
that the self-interaction tries to rearrange the free vacuum which was imposed 
as an initial state. It is possible to prevent this rearrangement if we choose 
constant C to be: 

This constant cancels exactly the self-interaction term and the initial free 
vacuum becomes a stationary solution to ( |103j - |106| ) and stays unchanged as 
long as t < to- Note that C is defined by the initial vacuum and stays fixed 
even when we switch on an external potential. 

The only problem is that the integrals defining C and self-interaction 
diverge and one has to introduce a finite cutoff A to make the subtraction 
meanigful. We work with an isotropic system so the angular integration is 
not a problem: 

//*oo /*A 
d 3 p -> 4vr / d\p\p 2 ^A-n d\p\p 2 . (110) 
Jo Jo 

The last substitution is just the desired regularization. The presence of di- 
vergences is not a strange here. In fact the constant C may be traced back 
to the "mass renormalization" counterterm in the Lagrangian. (There may 
be some confusion related to the fact that "mass" appears not only in U 
but also in other parts of our expressions. This is due to Feshbach-Villars 
transformation. The mass term plays an important role in this transforma- 
tion but we are free to use only a part of it - the rest may be attached to 
a scalar potential. Note, however, that our wave functions and subsequently 
the Wigner function are normalized with respect to the mass m i.e., the one 
that participates in the Feshbach-Villars transformation. This must be kept 
in mind when interpreting the results.) 

We solved evolution equations ( |104| )-( |rO"S| ) numerically for a finite volume 
of phase space < \p\ < A. This automatically introduces the cutoff. We 
have chosen it to be rather large so that we did not observe its effects. 

The non-trivial evolution starts at t — to when we switch on the — M 2 
term. It was chosen large enough to override the initial positive m 2 and 
the system rolls down towards the new ground state. The results are shown 
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in Fig. [I]. The maximum at low momenta indicates the development of 
long-range correlation as shown by Eq. (|101|) . (It is not exactly a delta 
function because we are dealing with finite times.) A more detailed analysis 
shows that the 'classical component' of the field reaches the value given by 
(|100 ) for large times when the oscillations around the new minimum die 
out. Both the initial decay and damping of oscillations are purely quantum 
processes (the latter results from particle creation by a time-dependent scalar 
potential). We should stress that this numerical analysis was intended only 
as an illustration and there are still some questions left. One of them is 
the stability of the procedure described here although we observe the good 
numerical stability of our results. 

The methods developed here allow for the study of far more complicated 
transitions in presence of abelian gauge fields, spatial inhomogeneities etc. 

Our formalism does not include the temperature effects which are of great 
importance for the phase transitions (e.g., in cosmology). This is an obvious 
shortcoming. We think that thermal effects can be incorporated at moderate 
expense. 



5 Summary 

Starting from the definition of the Wigner transform as a realization of the 
Heisenberg algebra in phase space, we have developed a calculus to trans- 
late operator equations into phase-space equations that have an expansion 
in powers of h and clearly exhibit the classical limit. Applying this formal- 
ism to the symmetrized correlation function of the Klein-Gordon field, we 
have derived the equation of motion for the equal-time Wigner transform of 
the Klein-Gordon field and given a meaning to its components. The result- 
ing set of equations has been shown to possess an interpretation in terms 
of a relativistic gas of negatively and positively charged particles. Interfer- 
ence terms gave rise to quantum corrections and, especially, Schwinger pair 
creation. This formulation of the Klein-Gordon field is especially suited to 
semiclassical nonperturbative calculations as is demonstrated in the case of 
spontaneous symmetry breaking. 
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Figures 
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Figure 1: The decay of the false vacuum. Initial mass m — 1, A = 1.0. The 
sign of the mass squared term is reversed at t = t = by switching on 
M 2 = 2. (Units defined byc = h=l.) 
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